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Abstract 

Gauge mediated supersymmetry breaking scenarios with an ultra-light gravitino 
of mass 7713/2 = 1-10 eV are very interesting, since there is no cosmological gravitino 
problem. We propose a new experimental determination of the gravitino mass for 
such an ultra-light gravitino, by measuring a branching ratio of two decay modes of 
sleptons. 



1 Introduction 



Gauge mediation of supersymmetry (SUSY) breaking in a hidden to the SUSY standard 
model (SSM) sector [1] is very attractive, since it provides us with a natural solution 
to the problems of flavor- changing neutral current and CP-violation in the SSM. The 
gauge mediation predicts a light gravitino (G3/2) of mass m 3 / 2 = 1 eV — 10 GeV, which 
is most likely the lightest SUSY particle (LSP). In particular, an ultra-light gravitino of 
mass in a range of 1 — 10 eV is very interesting, since there is no cosmological gravitino 
problem [2]. For such an ultra-light gravitino, the determination of its mass will be of 
crucial importance at the LHC and ILC, since the small gravitino mass would predict a 



relatively low SUSY-breaking scale ^J~F = \J V3m 3/2 M P ~ 60 - 200 TeV, where M P ~ 
2.4 x 10 18 GeV is the reduced Planck scale. If it is the case, it might not be impossible 
to access the physics for the SUSY breaking in future experiments. 

Nevertheless, it is very difficult to determine the gravitino mass for the light gravitino. 
The gravitino mass may be determined by measuring the lifetime of the next lightest 
SUSY particle (NLSP), since the NLSP lifetime is proportional to F 2 oc m\ j2 (HHjtl For 
the ultra-light gravitino, however, the NLSP lifetime is very short and it will be very 
challenging to measure its decay length: For a slepton NLSP, the decay length is given 
by ct ~ 0.55 /im(m 3 /2/l eV) 2 (mNLSp/200 GeV)~ 5 , whereas the experimental resolution 
of the impact parameter is at least 0(10) fim or larger [6j [7J. 

In this letter we propose a possible independent experimental determination of the 
mass for such a light gravitino £3/2- We consider the case where the NLSP is the lighter 
stau (f), and the right-handed smuon (/i) and selectron (e) are slightly heavier than the 
stau (f) but lighter than the lightest neutralino x\i m f < m ji/e < m x?' wmcn is a typical 
mass spectrum in the gauge mediated SUSY breaking scenario (GMSB) as shown in the 
next section. The smuon/selectron then decays into two channels; one is £ f + t + £ 
and the other is £ — * £ + G3/2, where £ {£) denotes fx or e (// or e). The former process 
is fixed by the SSM interactions and SUSY-particle masses and the latter is fixed by the 
gravitino mass and the smuon/selectron mass. Therefore, we can determine the gravitino 



x If the gravitino mass is large enough, a measurement of the gravitino mass by kinematical recon- 
struction may also be possible, which serves as a test of supergravity [5]. However, this is also impossible 
for the ultra-light gravitino. 
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mass by measuring the ratio of these decay processes, provided that the masses of relevant 
SUSY particles are known. 



2 Gauge-mediation model 

We consider a simple gauge-mediation model, where a SUSY breaking field S couples to 
N pairs of messenger chiral superfields, if) and ip, which transform as 5 and 5* under the 
SU(5)g\jt- W = kipTpS. The S field develops a vacuum expectation value k(S) = M + 
8 2 F, where M is the messenger mass. We assume an ultra-light gravitino of mass 1 — 10 eV, 
and therefore it is natural to assume a direct gauge-mediation, that is F/k = y^m^^Mp. 
We choose k — 1 in the following discussion, for simplicity. 

With the above conditions the low-energy spectrum of the SUSY particles including 
the gravitino mass are determined by 5 parameters, F, M, N, tan/3, and sgn(/i) = ±1. 
We consider iV < 4 to maintain the perturbative unification of the standard model (SM) 
gauge couplings. 

Let us now discuss the spectrum of the SSM SUSY particles. In the GMSB models 
gaugino masses are generated from loop diagrams of the messengers pp. At the one-loop 
level, gaugino masses are given by 

M a = ^Ag(x) (a = 1,2,3), (1) 

where A = F/M, x = F/M 2 , a t = 5a EM /(3 cos 2 9 W ) and 

g(x) = -^[(1 + x) log(l + x) + {l-x) log(l - x)}. (2) 

Scalar masses, at the two loop level, are written as 

^=2AA 2 E(^) 2 C a (z)/(x), (3) 

where C a {i) are Casimir invariants for the particle (pi (C\(i) = 3Y 2 /5) and 

(x -> -x). (4) 



log(l + x) - 2Li 2 (x/[l + x\) + -Li a (2x/[1 + x\) 



x z 

For x < 1, both f(x) and g(x) are 0(1). Detailed behaviors of these functions are given 
in p. F < M 2 should be satisfied, since otherwise messenger scalars are tachyonic and 
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Figure 1: m^-'s dependence on A for tan/3 = 
10. The lines with a short range are for 
m 3/2 = 1 eV, and the others for m 3 / 2 = 
10 eV. 
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Figure 2: m^o/m^ dependence on A for 
tan/? = 10. These ratios are almost inde- 
pendent on tan /3 



the SM gauge symmetry is spontaneously broken. We consider the case of F < 0.8M 2 in 
the present analysis. The masses in Eq.flTJ and Eq.([3]) are given at the messenger scale. To 
calculate the decay rate of the scalar leptons, we need to evaluate the on-shell masses and 
mixings for the SUSY particles. We should evolve the renormalization group equations 
down to the electroweak scale and calculate the neutralino mixing and scalar's left-right 
mixing. To evaluate those values, we have used the program S0FTSUSY 2.0.11 [8]. In our 
computation, we have chosen sgn(/x) = +1. 

For our purpose, m f , m~ t and m^o are the most important parameters. In Fig{TJ m~ t are 
shown as a function of A. We see that m^ depend linearly on A. We have checked that m^ 
are almost independent of tan/3. Figj2] shows the ratios m^o/m^. Those are also almost 
independent of tan/3, and change mainly through 0(1) functions, f(x) and g(x) in Eq.([T]) 
and Eq.flSj). For N — 1, the mass spectrum cannot satisfy our assumption (m f < m~ t < 
m^o). Therefore we consider the case > 2 in the present analysis. Lastly, Fig 13] shows 
the contour plots of Am = m^ — mf. We see that Am is a simply increasing function of 
tan/3. For Am < m T , the decay £^£-\-r-\-?is kinematically forbidden, and effectively 
all the I decay via I — > i + G3/2 for the ultra-light gravitino. Thus, we concentrate on the 
case of Am > m T , which corresponds to tan /3 ^ 5. 

We should also note here on the constraints from the lower bound of the Higgs mass. In 
large parameter regions of our interest, the so-called decoupling limit is satisfied. Namely, 
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Figure 3: The contour plots of the mass difference between m? and rug, Ara = m < - — m?, 
on the (tan/3, A) plane. The contour parameters are 0.5, 1, 2, 5, 10, 20 GeV from left to 
right, (a) is for iV = 2 and m 3 / 2 = 1 eV, (b) N = 4 and m 3 / 2 = 1 eV, (c) N = 2 and 
m 3/2 = 10 eV, (d) iV = 4 and m 3/2 = 10 eV. 



CP-odd Higgs A is heavy: m^o ^ 200 GeV. Thus, the lighter CP-even neutral Higgs hP 
is an almost SM-like Higgs particle. We find that some parameter regions are excluded 
by the experimental lower bound m/jO ^114 GeV. For example, all parameter regions are 
excluded in the case of N = 2 and m 3 / 2 = leV as shown in Fig 16]- (a). 

3 /2/e decay 

In our model fl/e can decay in two channels, as shown in Fig H] and Figj5j In the following, 
we neglect the masses of the electron and the muon. The decay rate T(£ — > i + f + r) is 
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Figure 4: A diagram for 
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Figure 5: A diagram for I — > £ + G3/2 



calculated from a Feynman diagram in FigJU where £ = fi or e and £ = /i or e. Detailed 
calculations are given in [9]. For the parameter regions of our interest, the decay rate is 
approximately written as 



T{£^£ + f + r) 



2 / a EM \ 2 Mm\ 5 



157T Vcos 2 6 



(4° 



m~ t 



44eV( m i V 1 + 

.100 GeW UOGeVy (r 2 -l) 2 ' 



(5) 
(6) 



where r^o = m^o / m~ t Here we have assumed Am = m^ — mf 3> m T . When Am becomes 
comparable to m r , T(£ — > £ + f + r) becomes even smaller than Eq.Q. On the other 
hand, T(£ —> £ + G3/2) is calculated from a diagram in Figj5l For a light gravitino, the 
decay rate is given by 

,5 

(7) 



T(£ 



G3/2) 



m- n 



487T m 2 /2 Mj, 



fm 3/2 



(8) 



We see that the decay rate T(£ 
other hand, the decay rate T(£ — > £ 



°' 011 eV {Tmtiv) {lev 

£ + G3/2) becomes larger for smaller m 3/ / 2 . On the 
f + r) is suppressed for smaller Am. Therefore, in 
a certain parameter region, these two decay rates become comparable to each other. We 
have found parameter regions of (N, tan j3, A) in which the ratios, T(£ — ► £ + G3/2) /r(£ — > 
£ + f + r), are of 0(1) (from 10 -2 to 10 2 ) as shown in Figj6l From Eq.Q, the gravitino 
mass is written as 

'T(£ £ + f + t)\ 1 



ml 



m 



3/2 



48ttMJ, v r(£ 



G 



3/2, 



r(£ 



r + r 



(9) 



If masses for relevant SUSY particles are known, we can calculate the decay rate T(£ 
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Figure 6: The contour plots of the ratio of the decay rates, T(£ — > G 3 / 2 £)/T(£ 
the (mj-, Am) plane. The contour parameters are 10 3 , 10 2 , 10, 1, 10 -1 , 10~ 2 
is for iV = 2 and m 3/2 = 1 eV, (b) N = 4 and m 3/2 = 1 eV, (c) N = 2 and m 3/2 = 10 
eV, (d) N = 4 and m 3 / 2 = 10 eV. Below the thick solid line, the assumption m^o > 
is not satisfied. With the constraint of m h o > 114GeV, all the region of (a) is excluded, 
and (b), (c) and (d) are partially excluded. 

e + f + r) ^ Therefore, we can determine the gravitino mass m 3 / 2 from Eq.Q by 
combining the measured branching ratio, T{£ — > £ + G 3 / 2 )/T(£ — ► £ + f + r), and the 
theoretical estimate of F(£ — > £ + f + r). 



2 For our parameter regions, 0.6 ^ T(l —> I + f + r)/T no m ixing(^ — * I + f + t) ^ 1, where r no mixing is 
calculated without the mixing parameters. Thus, we can determine TO3/2 with an accuracy of factor 20% 
without knowing the mixing parameters. 
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4 Summary 

Once SUSY particles are discovered at the LHC/ILC and their masses are measured, 
one of the most important tasks in the next step would be the determination of SUSY 
breaking scale, or the gravitino mass. In this letter we have proposed a new experimental 
method to determine the gravitino mass, which works well for an ultra-light gravitino of 
mass m 3/ / 2 = 1-10 eV. We have found that the branching ratio of T{£ — > G3/2 + £)/T(£ — > 
f + r + £) becomes 10~ 2 — 10 2 in a broad parameter region for the gauge-mediated SUSY 
breaking scenario with the gravitino mass m 3 / 2 = 1-10 eV. We can estimate the gravitino 
mass with the measurement of this branching ratio, provided that the masses of the other 
relevant particles are known. The largest error of the gravitino mass may come from the 
error of Am in our method. To determine the gravitino mass with an accuracy of factor 2, 
we need to know Am with an accuracy of factor 30% (cf. Eq.Q). This may be accessible 
by precise measurements of the slepton masses and/or by studying the decay products of 
the process £ —* r + r + I. 

Remarkably, we have also found that, the dominant decay modes of the smuon and 
the selectron are £ — * £ + G3/2 in a large parameter region. These processes produce an 
energetic muon or electron together with missing energy, which may become clear signals 
for an ultra-light gravitino scenario. A study of possible backgrounds in determining the 
ultra- light gravitino mass at the LHC/ILC will be given elsewhere |10j . 
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